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Abstract. Choice models for populations of agents on graphs are studied in terms of statistical thermo- 
dynamics. Equations of state are derived and discussed for different connectivity schemes, utility approxi- 
mations, and temperature and volume regimes. Analogies to ideal classical and quantum gases are found 
and features specific for network systems are discussed. 
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1 Introduction 

Ensembles of active, ehoice-making individuals represent 
a natural field of extension for statistical physics. This 
broadening of theory's scope is interesting for a physicist 
because of a new quality brought to the study by auton- 
omy of system's constituents in making decisions. In ad- 
dition, it is not a priori clear how to formulate an analog 
of the classical, equilibrium Boltzmann-Gibbs thermody- 
namics for systems where the notion of physical energy is 
not defined. How probability measure on the state space 
is defined then? 

Choice models, being part of decision theory, are widely 
applied to economics, psychology and social sciences. For 
example, models used in demand analysis of transporta- 
tion and communication usually assume that demand rep- 
resents the result of decisions of individuals in the popula- 
tion. These decisions consist of choices made among finite 
sets of possibilities. Consider decision to be taken by po- 
tential passenger willing to travel between two nodes. This 
client has to choose among offers of carriers, accounting 
for many factors relevant for decision, as e.g. the timetable 
and how it relates to his needs, durations of travels, fares, 
numbers of stops and changes on the route, probability 
of delays, declared and expected quality of service and 
many other aspects, too numerous to itemize. In order 
to quantify such process of decision making, one incorpo- 
rates discrete choice models in hope of better understand- 
ing and predicting behaviour of such complex system as 
transportation network and thus obtaining hints for mar- 
keting and revenue management. 

Provided one knows how to formulate statistical ther- 
modynamics for populations of agents, our understand- 
ing and predicting power for such systems is enhanced, 
as compared to choice models. Having probability dis- 
tributions in equilibrium one calculates thermodynamic 



potentials and response functions, and their evolutions 
with intensive variables as temperature, fugacity or pres- 
sure. Moments of extensive variables, directly related to 
thermodynamic potentials, have clear economic value, e.g. 
market and route occupancies or fluctuations of numbers 
of clients on them. In addition, monitoring collective phe- 
nomena and phase transitions is feasible using thermody- 
namic formalism. Such tools are lacking in economy and 
social sciences. 

In our previous papers we demonstrated how the sys- 
tem consisting of a set of decision makers on a network can 
be treated in the framework of equilibrium statistical ther- 
modynamics and its generalization using Renyi entropies 
PE]. Applying thermodynamic formalism one tacitly as- 
sumes existence of stationary regime and states optimal 
with respect to entropy. Finding thermal equilibrium is 
therefore related to the class of optimization problems, 
in particular to optimal path distribution on the network. 
Problems of optimal path search in discrete systems, quite 
often met in condensed matter, lattice field theory, net- 
works of molecular reactions, telecommunication and en- 
ergy routing, were often treated non-thermodynamically 
(cf. e.g. refs 0). We believe that application of Metropolis 
algorithms to this class of problems could be also effective. 
Interesting for us, however, was not the optimization itself 
but studying global properties of the system in quasistatic 
approximation, assuming ergodic hypothesis. For this we 
found thermodynamic approach to be well suited. 

Building thermodynamic description, the central point 
was to find a probability measure on the event space. 
We studied analogy between the utility function, used to 
quantify decision likelihood, and the total energy func- 
tion describing probability of states in physics. We found 
that for a properly defined system and utility function this 
analogy can be made exact. In particular, this concerns 
the properties of additivity and extensivity of utility, and 
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entails incorporation of the Boltzmann-Gibbs state proba- 
biUty measure, as for physical systems with no long-range 
correlations. 

Building this analogy, however, we have not made it 
complete and never arrived to the equation of state for 
populations of agents on graphs. In this paper, after re- 
capitulation of definitions of the system and statistical 
ensembles in chapter |21 wc introduce another extensive 
variable for the system, the market and network volume, 
and the conjugate intensive variable being an analogue 
of the pressure, and discuss equation of state of the net- 
work, both in general and for specific network topologies. 
Analogies to behaviour of classical and quantum gases are 
highlighted, wherever found. 

2 The system and the statistical ensembles 

Our system S consists of a network represented by a 
set of directed graphs, and a set of agents willing to find a 
directed route from the origin node O to the destination 
node D. Each ordered pair (OD) and set of all Hamilton 
routes^ from O to D constitute, adopting economic termi- 
nology, the fc-th market and is represented by a directed 
graph Qk {k — 1, . . . , M, M denoting the total number of 
ordered (OD) pairs in the network). The whole network 
is the union of market graphs, G = tJi U . . . U Qm- Market 
graphs can overlap, e.g. two subgraphs Qi and Qj can have 
common nodes or edges. Such decomposition of the mar- 
ket graph Q induces decomposition of the system S into 
the sum of markets iS = 5i © . . . © Sm, where © stands for 
the sum of graphs and agents ascribed to corresponding 
graphs. Fig.nillustrates decomposition of an example sys- 
tem into markets. Each market {OD) can be further split 
into routes [0,Ci, . . . ,Cl, D] in the same way, e.g. for the 
network from Fig. ^ decomposition is the following: 



1,2) 


= [1,2] 


1,3) 


= [1,2, 3] ©[1,2, 4, 5, 3] 


1,4) 


[1,2, 4] ©[1,2, 3, 5, 4] 


1,5) 


= [1,2, 3, 5] ©[1,2, 4, 5] 


2,3) 


= [2, 3] ©[2, 4, 5, 3] 


2,4) 


= [2, 4] ©[2, 3, 5, 4] 


2,5) 


= [2,3,5]© [2,4,5] 


3,4) 


= [3, 2, 4] ©[3, 5, 4] 


3,5) 


= [3, 5] ©[3, 2, 4, 5] 


4,5) 


= [4, 5] ©[4, 2, 3, 5] 




plus all reversed (OD) (DO) pairs 



The state of the system is defined by choices of routes 
performed by all agents for all markets. The state space 

^ A Hamilton route is a path between two vertices of a graph 
that visits each vertex exactly once. Restriction to Hamilton 
routes is reasonable for applications to communication and 
transportation problems where loops should be avoided. Gen- 
eralization of the formalism to all conceivable paths on a graph 
is straightforward, however we did not try to find analytic so- 
lutions in this case. 




M n M 

111 1 

full (1,2) (1,3) (1,5) 

network 
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(2,3) (2,5) (3,4) 

Fig. 1. Decomposition of the communication network into 
markets, where only examples of representative markets are 
shown. The full list of markets and their decomposition into 
Hamilton routes is given in the text. 



consists of all possible distributions of all agents for all 
markets. 

We consider so called disaggregate model where agents 
decide independently of each other. This assumption is 
not restrictive for a wide class of applications. Notion of 
an agent can be often redefined and refer to groups of 
agents taking decisions as a whole, and only mutual inde- 
pendence of groups is relevant. Generally, the mechanism 
of group decision making can affect utility function. We 
do not pretend to discuss rather vast subject of group de- 
cisions and refer to one of recent reviews ^ . 

The utility Ui^ of the ifc-th agent's choice {ik = 1, • . . , Vk, 
Vk being the number of choice alternatives on the k-th 
market and k = 1, . . . , M the market index) is a finite, 
real random variable on the choice set and is defined ax- 
iomatically using the preference- indifference operator 
This operator relates any pair of alternatives from the 
choice set and is assumed to be reflexive and transitive, 
thus ordering the choice set linearly. The choice set is as- 
sumed to be finite which guarantees existence of the best 
alternative. These properties of the preference-indifference 
operator induce the linear order in the real-valued utili- 
ties and ensure existence of the largest utility for the most 
prefcred alternative. Utilities are therefore upper-bouned 
and disutilities, defined as negative utilities, are bounded 
from bottom. 

The concept of utility dates back to the XVIII-th and 
XlX-th century utilitarian's economy and its eminent rep- 
resentative Jeremy Bentham. In its further evolution it 
acquired quantitative and strict meaning, mainly due to 
Morgentern and von Neumann |5] . In order to use it in spe- 
cific applications of choice models, one has to postulate a 
lot from outside of the theory and to find an effective way 
to either derive it from any fundamental or general theory 
or to parametrize it and estimate from data. For exam- 
ple, efficient estimates of utility in marketing is done in 
the framework of conjoint methods 0. Derivation of the 
utility function from first principles represents major the- 
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oretical problem because of lack of fundamental theory 
of human behaviour. Normally, one has to rely on par- 
tially justified models using power series approximations 
for utility functions, often restricted to linear or quadratic 
terms (cf. ref. [S]), and finding empirically relevant vari- 
ables and estimators for parameters, if any data are avail- 
able at all. For a broader discussion we refer to e.g. ref. 

Apart from those problems, numerous conceptual dif- 
ficulties arise when applying utility functions to quantify 
the behaviour of the choice makers and to predict them. 
In particular, classical Bentham's approach of maximizing 
the overall utility, integrated over individuals, when taking 
market decisions, is often a subject of serious objections 

nni- 

For each individual, utility is reduced whenever choice 
is performed, similarly to the wave-packet reduction in 
the process of quantum-mechanical measurement. Conse- 
quently in this paper we denote utility as a random vari- 
able by U and its particular instance, e.g. after making a 
choice, by u. For disutilities we adopt notation U and u, 
respectively. 

The system consists of TV = X^fcli -^fc agents, where 
Nk stands for the number of agents on the k-th market. 
Individuals' disutifities Ui^ = —Ui^ contribute additively 

to the overall network's disutility U = X]i2=i ■ If mar- 
kets are assumed to be atomic subsystems and agents do 
not interact then additivity of the overall system's disu- 
tility is ensured by construction [l]. Indeed, decompose 
the system S into the sum of subsystems (markets) Sk, 
S = iSi © . . . © Sm- Then, by construction, additivity 
condition is fulfiled U = IJi + . . . + Um (for more de- 
tailed discussion of additivity cf. ref. CJ). Furthermore, 
for U = 0{N°'), a < 1, disutility is also extensive. 

In order to follow classical Boltzmann arguments lead- 
ing to the exponential probability V^U — u) ^ exp(— /Ju), 
besides additivity we also assume conservation of disutil- 
ity for the whole system. Validity of this assumption is 
traditionally matter of debate (cf. e.g. In economy 

one normally does not contest it for utility being function 
of money or any other conserved commodity. 

We further assume existence of the stationary regime 
and existence of the first moment of disutility, (U) < oo, 
at least in the limit of long times. This ensures existence 
of the equilibrum temperature T^q = l//?eg in that regime. 

We consider the canonical and the grand canonical en- 
sembles of systems (cf. ref. H). In the canonical ensemble 
the total number of agents N is fixed and the partition 
function can be written as 



M 



(1) 



k=l 



where Zl (/3) stands for the partition function for the fc-th 
market with one agent on it and Nk = 0, 1, . . . , oo. Eqn. 

assumes non- subjectivity of the utility, i.e. that utility 
of given route is the same for any decision maker. For the 
grand canonical ensemble, or random N, one introduce 
another Lagrange multiplier, called the chemical potential 



/i, and the grand partition function reads 

M 



k=l 



(2) 



where — Z\{(3), v — /3/i. Moments of extensive vari- 
ables and the entropy are given by the partition functions 
Z and S', and their derivatives over /3 and v. 

In order to make thermodynamic description complete, 
we define the volume Vfe of the fc-th market as the to- 
tal number of choice alternatives available for each agent 
on this market, and the total volume of the system as 
V = X^fcli ^fc- Such defined total volume is additive by 
construction. Using market volumes, partition functions 
(jll2(l can be written as 

In Z(/3; AT, F; M) = ^ iV,. In ^ e''^"^-^ (3) 



and 



lnS(/3,j/;Ar,F;A/) 



M 

E 

fc=i 



ln(l - e"" 



, (4) 



where Uj^^ stands for disutility of the jfc-th agent's choice 
on the fc-th market, TV = (iVi, . . . , TVa/), F = (Vi, . . . , Vm)- 
We define the pressure as 



for the canonical ensemble and 



(5) 



(6) 



for the grand canonical ensemble. Such defined pressure is 
naturally interpretable thermodynamically as an effect of 
the size of the choice set on the free utility F 



P 



(-) 



(7) 



Since volumes are discrete variables, differentiation over 
them has to be defined with care. We adopt an approach 
of the finite-difference calculus for functions defined on 
an equally-distant grid with a unit space. Approximations 
for the TO-th order derivatives operators can be expressed 
using the forward (backward)-difference operators A (V) 



gr, 



-A' 



-V" 



71 = 1 



(8) 



where for any test function f{x) {x E Zi), Af{x) = 
fix + 1) - fix) and S7fix) = fix) - fix - 1). Using the 
forward- or the backward-difference operators expansion 
depends on which one gives more accurate approximation 
to the derivative's values and this normally depends on 
the distance of x from the domain's edge. 
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3 Partition functions and thermodynamic 
variables 

3.1 The cumulant expansion 

Eqns (|5I6|) express the pressure p via /3, v and volume V, 
and lead directly to equations of state. Depending on the 
temperature regime, one can expand partition functions 
around any characteristic temperature f3o and fugacity j/q 
of the system. The (3q and vq are determined by the mean 
utility and number of agents. 

The In Z and In S are generating functions for the cu- 
mulants of U and N . The In Zl can be expanded as 

lnZi(/3,/3o,H) - f2{{Ul!{p^f3o,Vk))) (9) 

where the cumulants of disutility can be related to ordi- 
nary moments (cf. e.g. ^H), e.g. 

m) = iu') - {u)\ 

((c73)) = (t/3)-3(i72)(c7) + 2(c7)3, 

etc. (10) 

For the grand canonical ensamble, expansion of S'^. reads 

In ='fe(/?,/?o,i^, I'o, Vfe) = 

OO 



n,m— 

N]:^{l3 = l3o.v = vo,Vu))) X 
(/3o - /?)" - yoT 



(11) 



where the bivariate cumulants can be related to bivariate 
moments |14) . e.g. 

{{UN')) = {UN)-{U){N), 
{{U^N^)) = 2{U){Nf - 2{N){UN) - {U){N^) 

etc. (12) 



3.3 Temperature evolution of the pressure 

As found in rcf. 2 , using Renyi entropies, defined for any 
g G M ( J standing for mmibcr of states) 



I T^lnE.tiPl, 9^1 
' l-ELiPjlnPi, 9=1 



(14) 



one relates partition functions at different temperatures, 
with q being scaling parameter 

In S{qP, qu, V) - qln S(/3, ly, V) ^ (1 - q)I,{p, V\\h) 

and similarly for Z . 

Using eqn. lO, this leads to more general equation 

p(g/3, qv) = qpiP, ^) + ^(1 - q f^'^gy , (16) 

being not only equation of state but determining also the 
temperature evolution of the pressure. 



4 Applications 

4.1 Networks of increasing path lengths: analogues to 
quantum-mechanical systems 

Assuming disutility function depending only linearly or 
quadratically on the topological (OD) distance, for Hamil- 
ton networks containing routes of all lengths one observes 
direct analogies to quantum-mechanical systems: quantum 
harmonic oscillator and the ideal quantum gas. 

Networks containing routes of all lengths may consist 
of one market (one {OD) pair) or many markets (more 
than one (OD) pair). Such classes of networks are very 
broad, even excluding non-compact and tadpole graphs. 
Examples of one-market and multi-market networks are 
given in Figs. |21 and O respectively. Dependence of the 
market volume V on the number of nodes L is a/L (upper 
left in Fig. 121), linear for one-market networks (lower left in 
Fig. 121), and quadratic for multi-market networks (upper 
left in Fig. which makes these networks rather easy to 
simulate. 



3.2 The high-temperature regime 

Applying the cumulant expansion around (3 ~ 1/T 
for fixed ly and calculating the pressure (O we get 



^9^/1 
^ dV \3 



914 / 1 {Nk) d{U) 



dV\l3 Vk dVk 
dM 

fi,y,N dV ' 



/3=0 



/3 9Var(i7) 
2 



dVk 



1 /ainS 

dM 



(13) 



For M = 1 the first term in eqn H13() corresponds to the 
equation of state of the ideal classical gas: pV (i = (N). For 
fixed M also dM/dV = 0. Corrections to the ideal case 
depend on the network topology and the utility function. 





Fig. 2. Example network topologies with one market and route 
lengths increasing linearly in steps of one. 
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Fig. 3. Example topologies for multi-market networks of routes 
with lengths increasing linearly m steps of one. 



Assume agent's disutility proportional to the number n 
of route segments, Un = uo+nui (n = 0, 1, . . .), where uq^i 
are constants. For one agent, this model exhibits analogy 
to the one-dimensional quantum harmonic oscillator with 
the energy spectrum En = hiu{l/2 + n). For N agents, 
additivity of disutility implies Un = uq + ui J2iLi ^i, n = 
(ni, . . . ,71^), and is the same as the energy spectrum of 
the A^-dimensional oscillator = hjj{N/2 + ni + . . . + 
njv). The values of constants "So = 0, ui = 1 are assumed 
without loss of generality. 

For any oscillator-like topology, excluding n — Q routes 
corresponding to tadpole graphs, the statistical sum reads 



lim Z{(3,V) 

V—^OO 



lim 

V — >OQ 



1 



6/3-1 
N 



3/3-1 



(17) 



for /3 > 0. Using first 3 terms of the expansion @ one 
gets the equation of state 



N 



N N 

y ^ 12 



(18) 



The first term of eqn. H18|) corresponds to the ideal clas- 
sical gas and it dominates at high temperature. The sign 
of correction to the ideal case is positive or negative, for 
PV > 6 or PV < 6, respectively, for which we hardly see 
any analogy in gas dynamics. At the equilibrium temper- 
ature P = ln[U/ {U — N)] the pressure is equal to 



N 



U/{U- 

n(i - 



for V > 1. 



(19) 



Considering random number of agents N , at the equilib- 
rium temperature /3 = ln[(C/)/((C/) — {N))] and equilib- 
rium fugacity V = \n[{N)^ /({U) - (iV))(l + (iV))], the 
equation of state reads: 



{N) 



m/m-{Nw -I 



(20) 



Let us assume agent's disutility proportional to the 
square of the route length, Un = U2n^. For the topologies 
of Fig. 121 and for one agent this model corresponds to the 
ideal quantum gas in a one-dimensional box and for the 
A''-agent disutility Un ~ U2 '^iLi '^f i ^he analogy extends 
to the A^-dimensional box. Depending on the number of 
agents allowed on each market, ~ 0, 1, . . . , oo or A^ = 
0, 1 (Pauli principle) and taking U2 = 1, one reproduces 
the case of the Bose-Einstein or Fermi-Dirac statistics, 
respectively, for which the grand partition functions are 
(cf. e.g. ref. UHl) 



lnS(/3,i.) 



-EnHl 

E„ln(l + e'^e 



/3C7, 



bosons 
fermions. 



(21) 



The Fermi-Dirac case can be interpreted as treating the 
whole group of agents as a single decision maker acting on 
the market (A^ = 1) or remaining passive (A^ = 0). More 
general dependence of t/ on n may also correspond to non- 
ideal quantum systems or different boundary conditions, 
or topology different than in Fig. |3| 

Formulae for Renyi entropies of quantum gases are 
given in refs. [2] and can be almost directly used for the 
evolution equation (|16|l . Explicit summation over states 
has to be performed, depending on the network topology 
(cf. Fig. 13). 



4.2 Maximum-connectivity networks 




Fig. 4. Maximum- connectivity network. 



For the maximum- connectivity or complete network (cf. 
Fig. 01 all M = L{L — 1) markets are topologically iden- 
tical. Within linear utility approximation one finds 



()! 



-Puil 



and 



L-2 



y, = (L-2)!^l 



(22) 



(23) 



For L ^ 1 and ui = 1 one approximates V = MVk — 
eL^Ll and, keeping only leading terms in L, arrives to the 
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equation of state 



P 



N 
0V 



N 



yinL' 



0, 



(24) 



For high temperature, f3 < \nL, and lowly populated net- 
work, N <^V\nL, the equation of state of the ideal clasi- 
cal gas is recovered from eqn. (|24|l . Correction to the ideal 
case is always negative, as for attracting molecular forces 
and quantum Bose statistics. 



itself. Similar behaviour wc observe in a physical system 
where the pressure depends on the way the volume was 
changed. Moreover, the pressure of physical gas depends 
on its statistical properties, e.g. ideal fermions exhibit 
higher pressure than bosons in the same conditions. In our 
case network topology affects the properties of the system. 
Two ways of evolution of network topology leading to the 
same final state are shown in Fig. |31 



4.3 Hub-and-spoke networks 







Fig. 6. Hub-and-spoke network. 



Fig. 5. Two sequences of networks evolving towards triangle 
network by decreasing numbers of nodes and keeping their topo- 
logical identity: maximum- connectivity network (lower) and 
harmonic oscillator network (upper). 



In general, equation of state for high temperature can 
be derived from eqn. For fixed number of markets, i.e. 
the number of markets independent of size of the choice 
set, dM/dV = 0, the formula for ideal gas is recovered. For 
variable M we find the volume V = eL{L — l)r{L —1,1), 
where r{., .) stands for incomplete Euler gamma function 
|13) . and accounting only for terms at 1/(3 in expansion 
(O we get 



TV 



+ 



TV 



dr{L- 1,1) 



13V /3e(2L- l)r2(L- 1,1) 



dL 



(25) 



Second term in (|25|) represents the finite-L correction to 
the ideal-gas term N/(3V and it vanishes for L ^ oo, in 
agreement with H24|l . 

Interesting property of network dynamics can be il- 
lustrated when calculating pressure in the simplest non- 
trivial case of the complete L = 3 graph, being at the same 
time the oscillator-like network with two routes, V ~ 2. 
The pressure depends on the model we choose: keeping 
terms proportional to the 0-th and —1-st powers of /3 in 
eqs. (|18|) and (|24|l one obtains two different solutions 



P 




maximum- connect ivity 
harmonic osciUator. 



(26) 



This is understandable because the pressure depends on 
the change of utility with volume and not on the volume 



For hub-and-spoke networks, being often met in local 
airline transportation or resource brokerage in Grid 
computing ^B] , there are only two classes of topologically 
different markets: Mi = — 3L -\- 2 spoke-spoke and 
M2 = 2{L — 1) spoke-hub, as shown in Fig.|H| 

The partition function 



Ml 
ki = l 



-2l3uiNki 



n 



(27) 



where TVi(2) ^re agents' populations on both types of mar- 
kets, does not depend on the volume V = L{L — 1) and 
thus the pressure is zero. However, the grand partition 
function {ui = 1) 



E{f3,y) 



1 



1 



(28) 



does depend on the volume and leads to interesting equa- 
tion of state. The complete formula is very lengthy and 
wc quote here only its large-volume limit 



P - 



(TV) , (TV)(l-e^) 



(3V /3y3/2 ' 
reducing further in the high-temperature regime to 
(TV) (TV) 



P 



(3V y3/2 



(29) 



(30) 



The same result can be obtained directly from eqn. \1'6\ . 
For (TV) ^ 1/3/2 formula H3U|) describes ideal classical gas. 
Positivity of the pressure {V > 0^) is ensured in this ap- 
proximation because /J^ <c 1 and V > 1. 
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5 Discussion and final remarks 



Systems of choice makers on networks can be treated ther- 
modynamically, using utihty function in a way analogous 
to the energy function and defining the system volume 
as the number of choice alternatives. Working within the 
framework of equilibrium thermodynamics we found equa- 
tions of state and discussed them in detail for specific net- 
work topologies, utility functions and in different temper- 
ature regimes. 

For a single-market and multi-market networks with 
linear utilities and route lengths covering all integer num- 
bers with no degenaracy. close resemblance to the one- or 
multi-dimensional quantum harmonic oscillator was found. 
Equation of state, besides the ideal-gas terms, contains 
corrections with no clear analogy in gas dynamics. For a 
quadratic utility function the system exhibits similarity to 
the ideal boson or fcrmion gas, depending on the number 
of agents allowed on one market. 

In case of the maximum-connectivity network, apart 
from deriving equation of state, we also discussed its de- 
pendence on the evolution of the connectivity scheme with 
volume. This property we found to be a particularly inter- 
esting example of how the network structure determines 
global properties of the system. It may suggest a possi- 
bility of phase transitions driven by variation of network 
topology and, we believe, deserves further investigation. 

We had a closer look at the hub-and-spoke connectiv- 
ity scheme as being important for applications. Analytical 
solutions for such system can be found and they simplify 
in the large-volume limit where corrections to the ideal- 
gas terms are negative. They are exponential with 1//3 
quenching in inverse temperature and depend on the vol- 
ume as 

In future work, more attention is certainly called for 
elucidating the role of possible correlations between agents, 
correlations between markets and of the form of utility 
function for this class of complex systems. Any mecha- 
nism giving dependence between subsystems, either ex- 
plicit correlation of choice alternatives |17) or interaction 
term between utilities jT8) , or limited node or link transmi- 
tivity, may result with non-additivity of utility and even- 
tually need for using non-extensive thermodynamics jl9| . 
In particular, it is natural to expect a non-linear rein- 
forcement of utility for the system of two agents mak- 
ing the same choice. Investigation of correlations between 
markets, either coming from explicit interactions between 
agents or from network constraints, should proceed along 
two lines. First, for weak correlations between subsystems 
one stays in the framework of Boltzmann-Gibbs thermo- 
dynamics and accounts only for agent interactions. Sec- 
ond, in case of unavoidable non-extensivity, the whole for- 
malism has to be reworked using power-law probability 
distributions. It would be very interesting if any nontrivial 
modification of the topology-dependence of the equation 
of state is observed, as compared to found in this work. 
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